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Summary 

Consider a string of n positions, i.e. a discrete string of length n. Units 
of length k are placed at random on this string in such a way that they do not 
overlap, and as often as possible, i.e. until all spacings between neighboring 
units have length less than k. When centered and scaled by the 

resulting numbers of spacings of length 1, 2,..., /c — 1 have simultaneously 
a limiting normal distribution as n —> cx3. This is proved by the classical 
method of moments. 
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1. Introduction and main result 

Imagine a row of n hooks on which one hangs hats at random. Each hat 
needs k adjacent hooks. Successive hats are put on the hooks by a random 
device that selects the next set of k adjacent hooks to be covered by a hat 
with equal probabilities from the remaining sets of k adjacent hooks. This 
process continues until all sequences of adjacent free hooks have size less than 
k. 

Let Xnj be the resulting number of sequences of adjacent free hooks of 
exactly size j. Necessarily j ranges from 1 to /c — 1. In the formulation 
of the problem given in the Summary, is the number of spacings of 

' j ‘ 

length j. In this paper asymptotic normality of X^ = (^n,i, • • ■, Xn,k-i) 
is proved. 

Theorem 1.1 Let X^ be the vector of numbers of spacings of length j, 
j = 1,..., /c — 1. Then 

(1.1) Zn = [Xn — EXn)-^ X'{0jT,) , as n —cx), 

and 

(1.2) [EXn — nO) —>0 ,asn—s>cx), 

hold. The components Oj, j = 1,..., k — 1, of 9 are described in Lemma 3.1 
and A/'(0, E) denotes the centered normal distribution with covariance matrix 
E = and with aij as dehned in (3.13) through (3.15). 

An outline of our proof of this limit theorem is given in Section 2. The 
limit behavior of expectation vector and covariance matrix of X^ is stud¬ 
ied in Section 3. The case k = 2 has been treated along different lines by 
Runnenburg (1982), who also gives references to literature about this prob¬ 
lem. As Runnenburg (1982) mentions, the random device we study here is 
relevant to problems in chemistry. Mackenzie (1962) describes its application 
to complete adsorption of linear molecules into parallel troughs (valleys) on 
suitable crystal surfaces. Here each trough is represented by our row of hooks 
and our hats symbolize the linear molecules. Mackenzie (1962) has derived the 
asymptotic behavior of the mean EX^ by a similar technique as in our Lemma 
3.1. The strength of our approach is in Lemma 2.1, which enables us to derive 
the asymptotic normality from the convergence of the hrst two moments. Our 
results are compared to these results from literature in Section 4. 

A continuous version of this problem has been suggested by Renyi (1958). 
It is called the parking problem and it has been studied by Dvoretzky and 
Robbins (1964). Coffman, Flatto, and Jelenkovic (2000) study the number of 
vacant intervals in this parking problem when the cars arrive according to a 
Poisson process in time and place. They consider the asymptotic behavior of 
mean and variance of this number of intervals at hxed time as the available 
space (the length of the row of hooks in the discretized version) grows beyond 
bounds. 
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2. Heuristics and a key lemma 

Clearly, when placing the hrst hat there are n — k-\-l possibilities, which 
all have the same probability l/(n — /c + 1). If the hrst hat occnpies positions 
j + 1 throngh j + k then two rows of hooks are left, one of size j and the 
other one of size n — k—j. Let J be a random variable nniformly distribnted 
on {0,1,..., n — /c} and let the {k — l)-vector Xi have the same distribntion 
as Xi. We assnme that J, Wq, Wi,..., Wq, Wi,..., are independent. 
The above rehection shows (cf. Lemma 2.4 of Rnnnenbnrg (1982)) 

(2.1) C{Xn) = C(^Xj + Xn-k-j) J n = k, /c + 1,.... 


In terms of characteristic fnnctions this can be formnlated simpler: 


( 2 . 2 ) 


Ee^ 




1 


i—k 


n — k + 1 ^ 

j=o 




s G n = k,k + l,.... 

With s = and as in the theorem this becomes 


(2.3) 


EE 


n 


-l n — k 

j=0 




) 


—j)/ Zn-k-j ^ g*”- ^^^X{EXj+EXn-k-j —EX^) 


Here Zq is degenerate at 0. 

If we conld prove that for some hnite constant B 


(2.4) \\EXj + EXn-k-j — EXn\\ < R, J = 0,..., n — /c, n = /c, /c + 1, .. . , 

holds with II ■ II the Enclidean norm and that {Z^} has a limit distribntion 
with characteristic f un ction Y then we conld dednce from (2.3) 

1 

(2.5) Yit) — J — ut)du, t G 

0 

By Lemma A.l from the Appendix this integral eqnation implies that 
{Zn} is asymptotically normal, bnt the asymptotic covariance strnctnre is not 
determined in this way, not even if we wonld know the limit behavior of the 
covariance matrix of Z^- Therefore, we will not try to elaborate the details of 
this approach, bnt we will nse the classical method of moments instead, which 
does yield the asymptotic covariance matrix. 
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This method of moments is based on the theorem of Frechet and Shohat 
(1931) abont the celebrated moment convergence problem, which states that 
if all moments of converge to the corresponding moments of Z and if 

the moments of Z determine its distribntion nniqnely, then {Z^} converges to 
Z in distribntion; see 11.4.C, p.l87, of Loeve (1977) or Takacs (1991). Onr 
key-stone for the method of moments is the fact that the normal distribntion 
is indeed determined nniqnely by its moments. By considering (c^Znj and 
c^Z for all c we see that the mnltidimensional resnlt is obtained from the 
one-dimensional resnlt. 


In Section 3 we shall show the validity of (2.4). Let Z^ 

= — EXn)^n = l,2,...,Zo = 0 a.s. Note that it follows from 

(2.3) and (2.4) that for all c G and all m G N, once we have 

snpj E \c^Zj\^ ^ < cx), then we obtain 


E{JZr.)^ 


( 2 . 6 ) 


1 


n—k 


n — k + 


15:^ 


j=0 


n 


n — k — j 


n 


7 1 

^ ^n — k—j 


H- -1= {EXj -|- EXn — k—j ~ EXn) 



m — 1 


i=l 


m 


1 


n — k / . X i/2 


i J n — k + 1 


Ui 


• E \ c 




3=0 
m — i 


n — k — j 




+ 


n 




E Z, 


E 


n — k + 1 ^ V n 

3=0 


E\ c" Z, 


as n —cx), 
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and similarly 

ElJZr, 


< 


1 


n — k 


n — k + 




j=0 


^'c^ZA + 


I^ ^n—k—j \ 


n 


(2.7) 


m — 1 

E 

i=l 


m 

i J n 


—— y.\- 

- /c + 1 ^ \n 

j=0 


+ \(^ {EXj + EXn-k-j — EXn) 


E c Z^—k—j 


+ 


n 

— k / . \ ml2 

J ^ 


E\c^Zj\ 


E 


n — k + 1 ^ 

j=o 


n 


E\JZ, 


These relations between the (absolnte) moments of enable ns to derive 

the asymptotic behavior of all moments of order m > 2 from that of the 
moments of orders 1 and 2, by the following key resnlt. 

Lemma 2.1 Let {anjJ^Lo ^ sequence of reals. If for some /c G N, a E R, 
and /? > 1 


( 2 . 8 ) 


ttn — CK + 


n — k + 


-lE. 




j=0 


n 


holds, then {a^} converges and 

(2.9) 


lim an ^ ot 

n^oo P “ 1 


Uj + o(l), as n —cx), 


/3 + 1 


holds. Moreover, if the elements of a possibly other sequence 

are nonnegative and bounded from above by the right-hand side of (2.8) for 

some /c G N, a G R, and /3 > 1, then snp^ is finite. 


Proof. If lim^^tx: 
i ^ oo, snch that 


\an\ — oo, then there exists a seqnence {n^}, Ui ^ oo as 


la^J = snp{|a/j| : h < nj, 


a. 


n-i I 


By (2.8) this yields 


1 = 


O'rii 

< 

O^rii 



oo. 


i-k 


< 


m-k + l ^ 

J=0 

r, rii — k 

-^y 

+ o(l), 


/ ■ \ 0 1 

( J ] 

ttj 

\nij 

l^n-i 1 


rii 


+ o(l) 


+ o(l) = 2 / x^dx + o(l) 


/3 + 1 
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which is a contradiction in view of (3 > 1. Conseqnently, we have 


( 2 . 10 ) 


lim \an\ < oo, and snp \an\ = a < oo. 


Incidentally, this proves the second statement of the Lemma. Let be a 

seqnence of integers tending slowly to inhnity; more precisely, oo and 

rrin/n —0 as n ^ oo. By (2.8) and (2.10) we obtain 


lim ttn <a + lim 

n^oo n^oo I n — /c + 1 


n-k / ■\ 3 


J=mr, 


Y. i snp ah + 2a 


nir 


nJ h> 


n — k -\-l 


( 2 . 11 ) 


n^oo 1^ J 

0 

2 

=a + —— hm a^. 

P + 1 n^oo 


snp ah 

h^m-n 


Conseqnently, we have 

( 2 . 12 ) 

Similarly, we obtain 


lim a„ < a 


/3 + 1 


lim a^, > CK + lim 
2 


_ 

n —^cxD n —^oo 


— Q) “t” 


/? + 1 


/3 + 1 h>m 

lim an, 


inf ah — 2a 


nir 


n — k + 1 


and hence 
(2.13) 


/3 +1 


lim > a ^ 

n —^oo P -L 


QED 


From (2.6) and Lemma 2.1 with (3 — mj^ > 1 
converge once the hrst two (absolnte) moments do 


it follows that all moments 
and more explicitly that 


(2.14) 


lim E 

n —^oo 



m 


m + 2 
m — 2 



^)(m-z)/2^^ 


■ hm E 

n —^oo 



lim E 

n —^oo 



m — i 


holds. Note that existence (and bonndedness in n) of the odd moments is ver- 
ihed here by application of the second part of Lemma 2.1 to the corresponding 
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absolute moments satisfying (2.7). Of course, the first moment vanishes. As¬ 
sume 


(2.16) 


lim E{c^Zn)‘^ = a‘^ < oo 


and note that this implies the existence and boundedness in n of the hrst 
absolute moments. Denote lim Eic^Zn) by Urn for m > 2. From (2.14) 
and by induction on m we arrive at 


(2.16) 


0 


odd 


fXrn = ■^ if m is 

2~'^ /mla'^ j{mj 2)1 even, 


which are the moments of a normal distribution with mean 0 and variance 
cT^. This yields asymptotic normality of c^Z^ by the method of moments 
as discussed in Section 1, since the normal distribution is determined by its 
moments; see e.g. Feller (1971), Examples VIII.6(b) and Vlll.l(e), pages 269 
and 251, 252. Consequently, the proof of (1.1) is complete once we have proved 
(2.4) and (2.15) with = c^Sc. This will be done in the next section. 


3. Asymptotics for the first two moments 

Asymptotically EX^ behaves as follows. 
Lemma 3.1 With the notation 


ek{y) =exp{2[y-SyV2 + --- + 2 /'' V(^-l)]} 

and 

1 

9j = 2{ek{l)) ^ j{l -y)y^ek{y)dy, j = 

0 

there exists for every R> 1 a constant Bn with 

(3.1) \EXnj - (n +k)9j\ < j = l,...,k-l, n = l,2,.... 

Note that (3.1) implies (2.4) and (1.2), and enables us to conclude from 

(1.1) that also 

(3.2) n-^/^(xn-n9]^J\f(0,i:) 


holds. Therefore, the proof of Theorem 1.1 is complete once (1.1) has been 
shown. 
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Proof of Lemma 3.1 Fix c G ^ and let 7 ^ = Ec^Xn- By (2.1) or (2.2) 


(3.3) 


and hence 


In 


1 

n — k + 1 


n — k 

5 ^( 7 j + In-k-j) 
j=0 


2 

n — k + 1 


n — k 


j=0 


n = k,k + 1,.. 
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(3.4) (n -/c + l) 7 n = (n -/c) 7 n_i + 27 „_fc, n =/c + 1,/c + 2,..., 


hold. Using generating fnnctions similar to Overdijk (1981) (see also Mackenzie 
(1962)) we dehne 


(3.5) 


Giz) = 


00 

n=fc+l 


for \z\ < 1 (note 7 „ = 0{n)). Mnltiplying (3.4) by ^ and snmming over 
n = /c + 1, /c + 2,... we obtain 

(1 — z)G'{z) = 2z^ ^G{z) + 2[7i2 + ... + 'jkZ^] + 'jkZ. 


In fact, 7 a; = 0 since = 0 a.s., and 'yi = Ci for i — 1,... ,k — 1. Solving 
this differential eqnation nnder the side condition G(0) = 0 we obtain 

G{z) = (1 - z)~‘^^{z), 

Z 

i^iz) = 2{ek{z))~^ j [ciy + ... + Ck-iy’^~^] (1 - y)ek{y)dy. 

0 

Note that U( 2 ), 2 G C, is an entire fnnction and write Giiz) — EfL^aiz'^. 
Then (3.6) yields 


00 n—/c+1 

(3.7) G{z)= ^ (n-ifc + 2-f)aiz”-^+^ 

n=fc+l i=2 

and hence by comparison to (3.5), for n = /c + 1, /c + 2,..., 

n—fc+l 

7n= ^(n-/c + 2- i)ai 

i=2 

00 CXD 

= ^^(n —/c + 2 — f)aj — {n — k + 2 — i)ai 

i=2 z=n—/c-|-3 

= [n - k + 2 )^( 1 ) - ip'{l) + Rn 


(3.8) 



with 


OO OO 

(3.9) |i?n| = I ^ {i - {n - k + 2))ai\ < ^ i\ai\. 

i=n—A:+3 i=n—A:+3 

Since '0(2), 2 G C, is an entire fnnction, ztjj^z) = entire 

as well and by Canchy’s root test for the radins of convergence of analytic 
fnnctions we obtain 


(3.10) lim = 0. 

i —^00 

Conseqnently, for arbitrary R > 1 and n snfficiently large we have 


(3,11) 


OO OO 

i\ai\ < R~^ 

i=n i=n 


R-^ 

l-R-^' 


Finally, note 

(3.12) {n — k + 2)'0(1) — 0^(1) — {n + /c)0(l) = (n + k)c^9, 

which by (3.8),(3.9), and (3.11) implies the lemma. QED 

For onr stndy of the (asymptotic) covariance strnctnre of we need 
some extra notation. Let 

Z 

(3.13) Gi{z) = 2{1-z)~‘^{ek{z)) ^ Jy\l - y)ek{y)dy 

0 

be the generating fnnction as in (3.5) throngh (3.6) of the means of Xn^i , i — 
1,..., /c — 1, -n- = /c + 1, /c + 2,..., and dehne 

Hij{z) =(1 - z)z^l[i=j] + (1 - 2 )^( 2 * + 2 ^“^Gi( 2 ))( 2 ^ + z'"~^Gj{z)) 

- - 2)-2 1^3 + (4k - 5)(1 - 2 ) + 2{k - 1)2(1 - z)2 

- 2/c2(1 - z)^ - (2 + {4k - 3)(1 - 2 ) + {2k - 1)^(1 - z)^ 

-4e{l-zf)z'^y i,j^l,...,k-l. 

We will show that 

1 

(3.15) aij = 2{ek{l)) ^ J Hij{y)ek{y)dy 

0 

is well dehned and even that 
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Lemma 3.2 For i, j — 1,..., k — 1, and for every R > 1, 

(3.16) Cov{Xn^i, Xnj) — {n + k)aij + 0{R~^), as n ^ oo, 
holds and consequently E of Theorem 1.1 satishes 

(3.17) E=(<Ty)‘;2j. 


Proof. Fix i,j G {1,..., /c — 1} and let = EXn^iXnj — {n + k)‘^9i9j for 
n = 0,1,.... By (2.1) or (2.2) we obtain 

2 n—k 

(3.18) Cn — ^ ^ ^ '^^{ch + {h+k)^9i9j+EXh^iEXn-k-h,j} — {'n^Fk)^9i9j 

and hence 

(3.19) {n — k + l)cn = {n — k)cn-i+2cn-k+‘2Sn,i,jj n = k + 1, k + 2,..., 
with 

n—k — l 

^n,i,j ~ ^ ^ {^h,i9j F 9iehJ “1“ ^h,i{,^n—k — h,j ^n—fc —/i —) } 

h=0 

F k9ien—k,j: 

^h,i EXfi^i (h k')9^, i 1,..., /c 1, /i 0,1,.... 

Tedions compntation shows 

OO CX) CX) 

n=fc+l h=0 n=k+h-\-l 

I h n — k — h h y'n—k — h — l '[ 

I Cfi^iZ t.fi—k — h,jZ tYi—k — h — l,jZ Zj 

OO 

F k9i ^ euqz'^ 

h=i 

=(1 - z){F F z'^-^G,{z)){z^ F z’^-^Gj{z)) 

F 9,9,{e - (1 - ^)-3(l + {k-l){l- ^))2}. 

By an argnment similar to that to derive (3.6) throngh (3.8) from (3.4) we 
obtain from (3.19) 

CXD 

(3.20) Gij{z) = ^ = {I-z)~‘^%l}ij{z) 

n=fc+l 
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with 


(3.21) 'il)ij{z) = 2[ek{z)) ^ Hij{y)ek{y)dy = ^bhz’^, 


Cn^ {n-k + 2)'0ij(l) - + Rr. 

= (n + k)aij + Rn 


0 h=3 

and consequently, as in (3.6) through (3.9) and (3.12), 

(3.22) 

with '0ij(l) = aij given in (3.15) and with 

OO 

(3.23) \Rr,\ < 

h=n — k-\-3 

To complete the argument along the lines of the end of the proof of Lemma 
3.1, i.e. along the lines of (3.10) we have to show that '0ij(2), 2 ; G C, is 
entire. Indeed, it may be verihed that, as 2 ; —> 1, the function Gi{z), z E C, 
behaves like 6*^(1 — 2 ;)“^(l + 2(/c — 1)(1 — 2 :)) + (!l(l) and hence Hij{z) like 0(1). 
Consequently, Hij{z), z G C, is entire, since (1 — z)‘^Gi{z) is. Together with 
Lemma 3.1 this implies that (3.22) yields (3.16). QED 

Note that by the argument of Section 2, (3.16) implies that all moments 
of c^Zn with Zn as in (1.1) converge to those of a normal distribution. 
Consequently, Theorem 1.1 holds with E as in (3.17). 

4. Comparison to literature 

For k = 2 and i — j = 1, Lemma 3.1 holds with 

(4.1) 6'i=2e“^/' {l-y)ye^ydy ^ e~‘^. 

Jo 

Furthermore, lengthy computations show 

(4.2) Gi(2) = (1-2)-V'-1, 


(4.3) 

and hnally 

(4.4) 


Hii{z) ^z(l ~z) + (l- 2(1 - z)-i + (1 - z)-2)e-^" 

- e“^((l - z)-^ + 2(1 - z)-^ + 1 + 29(1 - z) 

- 49(1 - zf + 16(1 - zf), 


an = 2e 


-2 


Hii(y)e^^dy = 4e 


-4 
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These results agree with the expressions for the asymptotic mean and variance 
obtained in (3.20) and (3.27) respectively, of Runnenburg (1982). Mackenzie 
(1962) has studied asymptotic mean and variance of the total vacant length 

k-i 

(4.5) = 

i=l 

Our Theorem 1.1 yields asymptotic normality of the standardized Vn- More¬ 
over, Lemma 3.1 and 

(4.6) 

fc-i j 

efc(l) = 2 (1 - y)iy + 2y‘^ + ... + {k - l)y'"~^)ekiy)dy 

j=i i 

1 

= 2 J {{ky -{k- l)y^)il + y + ... + y^~‘^) - (k - l)y)ekiy)dy 
0 

= [{ky -{k- l)y^)ek{y)]l -k [ ek{y)dy 

Jo 

1 

= efc(l) - k j ek{y)dy 
0 

yield the existence for every R > 1 of a constant with 

1 

(4.7) \EVr, -{n + k){l- k{ek{l))-^ J ek{y)dy\< Br R-^, n = 1,2,... . 

0 

This result is in line with (A28) and (A25) of Mackenzie (1962). An asymptotic 
expression with an exponentially small error as above can be obtained also for 
the variance of Vn via Lemma 3.2. 


Appendix 

Lemma A.l Any random variable on with characteristic function 'ip 
satisfying integral equation (2.5), i.e. 




1 

J^p{^/ut)^p{^/l — ut)du, t G 
0 


has a normal distribution. 
Proof. Fix a G R'^ and dehne 


(Al) 


Xa(s) 



e a)dw , s > b. 
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By the boundedness of ip this Laplace transform is well dehned and by domi¬ 
nated convergence we obtain 


(A.2) 


Xais) = ^Xa{s) 


e ^'^wip{^/w a)dw. 


Consequently, (2.5) yields 


(A3) 


~Xa('5) = / e / ip{y/uw a)ip{\/ (1 — u)w a)du dw 

Jo 

nW 

/ e~^'^ip{y/v a)ip{^/w — V a)dv dw 


'0 ^0 
/*oo /*oo 


'0 J V 

. 2 , 


^ — s(w—v 


^ip{^/w — V a) dwe ^'"ip{y/v a) dv 


=Xa(s) 

and hence there exists ba & C with 

Xa(s) = 


1 


ba + S 


, s > 0. 


By Laplace inversion (cf. e.g. VIL6.6, p.233, of Feller (1971)) this implies 

ip{y/w a) — ,w > Q. 


Consequently, a random d-vector Z with characteristic function ip satishes 

^ , s > 0, 

and hence Z is normal for any a G and Z is multivariate normal indeed. 
QED 
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